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Верификация модели пласта-коллектора 

по данным гидродинамических исследований 
скважин с использованием доверительных 
интервалов 


Представлена базовая процедура нелинейного регрессионного анализа. Использован байесовский вывод, 
при котором вся информация о параметрах пласта выражается в терминах вероятности, а неопределенность, 
связанная с оценками этих параметров, может определяться количественным образом. Исследованы 
проблемы, возникающие при использовании доверительных интервалов в рамках байесовского вывода 
для дискриминантного анализа различных моделей-кандидатов. 


Алгоритм нелинейной регрессии 


Цель нелинейного регрессионного анализа на основе метода наименьших 
квадратов (МНК) состоит в минимизации суммы квадратов разностей между изме- 
ренными и вычисленными значениями давления в рамках выбранной модели пласта: 


Е=>.(и-Р(и, м). (1) 


где ЕЁ — целевая функция; Г -— функция, описывающая модель пласта; 0 — неизвестные 
параметры пласта; х, — зависимая переменная (время); у, — независимая переменная 


(давление); и — число данных. 

Функция РЁ, а значит, и целевая функция Ё -— это обычно нелинейная функция 
неизвестных параметров пласта. Вследствие этого нужно многократно обновлять неизве- 
стные параметры до тех пор, пока целевую функцию нельзя будет сделать еще меньше. 
Метод Гаусса-Маквардта со штрафными функциями и специальной техникой интер- 
поляции и экстраполяции представляет собой модификацию метода Ньютона. Мате- 
матической базой метода Ньютона является теорема Тейлора. В методе Ньютона 
целевая функция аппроксимируется квадратичной моделью путем отбрасывания членов 
ряда Тейлора вблизи начального приближения набора неизвестных параметров: 


х -Еь +29, 6 „+2, 0, (6, 6, 9, . (2) 
где 
Ех). 
90, 2> (у Е(0 5 и ы 
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Е | 32| ое дР о Е 
00 00| — 00 00 ы Е(®', *)) 55 9 
та: 7 ке СЗ | 


Обозначим градиент целевой функции как © и матрицу Гессе как Н, которые 


о 


определяются как 


ОЕ 
=|о.|=| — |, 5 
Я 06 | (5) 
9 
0?Е 
Н=|Й, |=| ——| |. 6 
[Ин 00,09, | (9) 
В результате уравнение (2) может быть записано в матричной форме 
. 1 
Е Еф +(88) 9+5(58) Н(50), (7) 


где би=и-и’. 
Здесь Е” — аппроксимация по методу Ньютона значения Е в окрестности и’. 


Минимизация Е” требует, чтобы ее производная относительно би равнялась нулю в 
неподвижной точке: 


ОЕ" 
= 0. 8 
д(би) (8) 
Подстановка (7) в (8) приводит к следующему: 
би=-Н . (9) 


би, полученное из (9), является просто условием экстремума для Е”. Поэтому 
требуются итерации для нахождения решения для Ё. 
Новое решение определяется по формуле 
и=и° + би. (10) 
Итеративная процедура повторяется, пока не будет достигнута желаемая 
сходимость. Необходимо учесть, что когда би станет пренебрежимо мало, то значение 


Е” в неподвижной точке практически совпадает с Ё, то есть целевая функция до- 
стигает экстремума. С методом Ньютона связаны две очень важные проблемы, 
которые возникают в зависимости от природы функции, используемой для описания 
модели. Первая трудность касается членов, содержащих вторую производную, которые 
расположены на главной диагонали матрицы Гессе. В результате нет гарантии, что 
матрица будет положительно определенной, а значит, совсем необязательно, что на 
каждом шаге итерации новое решение будет стремиться к точке минимума. Вторая 
проблема связана с тем, что итерационный процесс может сходиться очень медленно 
или даже расходиться в случае, если матрица Гессе плохо обусловлена. Причинами 
могуг быть сильная корреляция, наблюдаемая между некоторыми параметрами, и не- 
чувствительность функции, описывающей модель, к определенным параметрам. Поэтому 
необходим ряд модификаций метода для преодоления этих трудностей, чтобы обеспе- 
чить быструю сходимость к точке минимума. В методе Гаусса члены, содержащие 
вторую производную, отбрасываются, как если бы они были равны нулю. Обычно 
такая аппроксимация является удовлетворительной из-за нулевого градиента целевой 
функции в точке минимума. Такая модификация делает матрицу Гессе положительно 
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определенной и гарантирует сходимость к точке минимума. Метод Маквардта полезен, 
когда матрица Гессе в методе Гаусса плохо обусловлена. Добавление константы к 
диагональным элементам этой матрицы улучшают ее обусловленность и препятствуют 
численной вырожденности матрицы Гессе. Специальная техника интерполяции и экстра- 
поляции изменяют длину шага для повышения скорости сходимости. Этот способ 
известен также как линейный поиск. Штрафные функции могут использоваться для 
дополнительного улучшения скорости сходимости путем ограничения области поиска. 

В итоге значение целевой функции в точке минимума вычисляется следующим 
образом: 


Е 


ЛЕ = 


„Е. +) $, -6,) 6) (10) 
00, |. 


и 


Это выражение эквивалентно следующей записи: 
"| ОЕ . 
Е(и)=Е(и) + -—| -(0,-6,) (12) 


Статистический вывод 


Со статистической точки зрения неявно предполагается, что функция, описыва- 
ющая модель, имеет следующую форму: 


у, =Е(и, х,)+ =, (13) 
где =, — случайная ошибка, характеризующаяся своей функцией плотности распреде- 


ления вероятностей. 

Оценка по МНК широко используется для решения обратных задач, т.к. по 
центральной предельной теореме сумма независимых случайных величин с конечной 
дисперсией стремится также быть распределенной по нормальному закону. Для нели- 
нейной модели вышесказанное уже, в общем-то, неверно. Однако считается, что для 
области в пространстве параметров достаточно близкой к оценкам нелинейную функцию, 
описывающую модель, можно аппроксимировать линейной функцией путем разложе- 
ния в ряд Тейлора первого порядка. Это допущение остается справедливым для 
большинства моделей пласта до тех пор, пока используется большой объем данных. В 


итоге для больших и даже в случае нелинейной регрессии оценка по МНК и может 
выступать в качестве несмещенной оценки и с минимальной дисперсией. 


Доверительные интервалы 


Непосредственное применение доверительных интервалов к результатам гидро- 
динамических исследований скважин (ГДИС) требует выполнения двух условий. Во- 
первых, ошибки, которые представляют собой разность между действительным значе- 
нием давления и его истинной величиной, должны быть независимы и нормально распре- 
делены относительно истинного изменения давления. Это предположение образует 
статистический базис МНК. Второе условие состоит в том, чтобы для области в про- 
странстве параметров достаточно близкой к их оценкам целевую функцию можно было 
аппроксимировать линейной формой путем разложения в ряд Тейлора первого порядка. 
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Условие возможности аппроксимации функции, описывающей модель, путем 
разложения ее в ряд Тейлора первого порядка приводит к следующему выражению: 


ОЕ _(©,-8,} (14) 


и. И=Ии 


Считается, что наблюдаемые данные замеров давления у, нормально распреде- 


лены относительно истинного значения Е(и, | с известной дисперсией: 


РиоБ( у Е (и, х,) = Роб (у; 


1 1 2 
и =" Е(и, х,)) } (15) 


В результате наблюдений и значений давления функция правдоподобия для 
параметров имеет вид 


и) = РиоБ(у, и)... РиоБ(у, 


н)- 
= ! 5 Е бы} = (16) 


1 2по 


тает вые [- (в. (ми) (®-5- (м-н) 


Киу,,..., у, = РиоБ(у,,..., У, 


где 
У Е(и м) 
В = : (17) 
У, -Е(и, х, 
и 
ОЕ ОГ 
|=), — [®] 
Т= : `В : : (18) 


И=и 


МНК эквивалентен максимуму функции правдоподобия, который имеет место 
тогда и только тогда, когда 


ВУ =0. (19) 


В результате 


Е(иу,...у,)= | 5 изв + (ии) зи) | (20) 


(М2по] 
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Матрица Гессе в методе Гаусса, поделенная на 2, определяется как 


Е НЕЕ 


Н= : 28. : : (21) 
[ОР || ОЕ п ОР ОЕ 


УУ=Н. (22) 

Если для параметров используется локально равномерное априорное распре- 

деление вероятностей (неинформативное априорное распределение вероятностей), 

тогда по теореме Байеса апостериорное распределение вероятностей параметров 
после и наблюдений: 


И=Ии 


Тогда 


ур ‚... У, |: РгГОБ уй оды 
а АН СЯ 
[Циу.... у, - РкоБ(и)4и [и...., у, аи 
где РгоБ(и) — локально равномерное априорное распределение вероятностей. 
По определению многомерного нормального распределения: 
О 
= у ехр ро (и и) н(и _ ") ам =1. (24) 
Поэтому уравнение принимает вид 

|н 1/2 1 АЯ ‚ 

Рио (Иу,,..., У, = - т ехр РЕЯ (м и) н("-*) (25) 


То есть параметры и образуют многомерное нормальное распределение отно- 
сительно и с матрицей ковариаций с”Н"'. Уравнение (25) количественно выра- 
жает неопределенность, связанную с оценками параметров. 

Когда дисперсия о” неизвестна, вышеописанные рассуждения требуют неболь- 


шого уточнения. о . можно получить на основе среднего квадрата ошибок 5 . р который 
вычисляется как 


а (26) 
п-т 
где 
5$ = У -Е(и,х,))'. (27) 


1=1 
В этом случае 5’ - это несмещенная оценка с”, и о’ имеет обратное гамма- 
распределение относительно 5” с п-т степенями свободы: 


2 и з 
Риов(с |5?) = Е | ый Е } (28) 


Поскольку ии о? — независимые случайные величины, то и не меняется, 
даже когда о? заменяется 52 Е 


где у=и-т. 
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В итоге апостериорное распределение вероятностей для и может быть получено 
путем исключения с” при интегрировании суммарной апостериорной плотности 
распределения вероятностей для и и с”: 


Роб (и| у»... у, )= [ Рио (и, в у, ао? = 
0 
Е [Реоь(щу,.-- р с”) Риов (о |.,..., у, ) о? = (29) 
0 


= [Рюов (му, а с”). Риов(в? |5? 49°. 
0 


После подстановки (25) и (28) в (29) получим 


гон" |1, (ия) н(и-#) | ° 


(Г(у2))"г(5/2)( ) у5 


12 


Рио (иу,,..., у,}= (30) 


Поэтому, когда с” неизвестна, параметры и образуют многомерное {-распреде- 


ление Стьюдента относительно и с матрицей ковариаций 5’Н" и п-т степенями 
свободы. 
Маргинальное распределение вероятностей параметра 9, определяется путем 


исключения 0, (1=/], 1=1,..., т) при интегрировании по пространству и: 


1 1 6 2 
Р’оБ(0|у,,...› у, == ехр Е (е,- й у (31) 
2поь. 0%, 
где о, — СКО, определяемое как 
в = бл (32) 
где й, — /-й диагональный элемент обратной матрицы Гессе, вычисляемой в точке 


^ 


И=и. 
По определению, 95 %-й доверительный интервал накрывает 95 % площади под 
кривой плотности распределения вероятностей, т.е. он представляет собой диапазон, 
доверительная вероятность попадания значений параметра внутрь которого составляет 
95 %. Обычно используются два типа доверительных интервалов: диапазон абсолютных 
значений и диапазон относительных значений. Относительные значения получаются 
путем деления абсолютных значений на величину оценки параметра. 


В случаях, когда дисперсия с” неизвестна, (1-а,)-100 %-й доверительный 


интервал для каждого параметра определяется из следующего неравенства: 


^ 


9, —о, о 39, <0;+5, “Па, (33) 


где (_„/› — табличное значение квантиля порядка 1 — 0/2 для {-распределения Стьюдента 


с и-— т степенями свободы. 
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В случаях, когда и-т> 30, значение 1_„», может быть заменено соответству- 
ющим значением для нормального распределения. Так, для @= 0,05 его значение 
будет равно 1,96. Тогда (33) примет вид 

0, —1,96-с°, =0,<0,+1,96-в... (34) 

(1 _ о) 100 % -й доверительный интервал для относительных значений каждого 


параметра определяется из следующего неравенства: 


О“ 0. бу ‘Па 
ее За р Е: (35) 
9; 9, 9; 


Аналогично, когда и-т>30, доверительный интервал для относительных 
значений каждого параметра можно представить в виде 
1.96.5, 0, 1,96-в.. 
< —<|1+ А. 
9, 9, 9, 
На практике используются значения доверительных интервалов, приведенные 
в табл.1 [1]. Они были получены эвристическим путем на основе реальных 
экспериментов по интерпретации промысловых и смоделированных данных ГДИС. 
Ключевой идеей является то, что если модель выбрана правильно и присутствует 
достаточное количество данных, то все параметры должны быть внутри этих 
допустимых пределов. 


(36) 


Таблица 1 — Приемлемые доверительные интервалы по [1] 


Диапазон Диапазон 
Параметры относительных абсолютных 
значений значений 
Проницаемость (А) +10 % 
Коэффициент накопления (С) +10 % 
Расстояние до границы (и, ) +10 % 
Длина трещины (х,) +10 % 
Коэффициент упругоемкости +20 % 
(©) 
Коэффициент пропускания ( /.) +20 % 
Скин-фактор (5) +1 
Первоначальное давление (Р) 0,005 МПа 


Продемонстрируем, как можно использовать доверительные интервалы для 
оценки корректности модели. Для этого были смоделированы данные ГДИС 
методом понижения уровня. Цель демонстрации — показать, как доверительные 
интервалы решают задачу, когда заранее известно, соответствует ли модель пласта 
данным или нет. 

В первом случае модель была выбрана правильно. Значения давления при ГДИС 
методом понижения уровня рассчитывались с использованием модели течения в 
бесконечном пласте, к которым затем прибавлялись случайные ошибки. Информация 
о пласте и насыщающем его флюиде представлена в табл. 2. Истинные значения 


параметров А=0,05 мкм”, 5 =10 и С=0,2 м’ /МПа . Генератор случайных чисел 
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формировал набор случайных ошибок, распределенных по нормальному закону с 
нулевым математическим ожиданием и дисперсией 2,5.10° МПа’. В зависимости от 


числа точек данных рассматривались следующие четыре случая: а) 51 точка данных, 
6) 61 точка данных, в) 71 точка данных и г) 81 точка данных. Применялась модель 
течения в бесконечном пласте с тремя параметрами (А, 5 и С). Соответствие модели 
данным иллюстрируется на рис. 1. 


Таблица 2 — Информация о пласте и насыщающем его флюиде 


Радиус ствола скважины (7. ) м 0,1 


Мощность пласта (й) м 5 
> 3 3 
Объемный коэффициент ( В. ) [м |" ) [м | 1 


Вязкость (и) Па.с 10° 
Пористость (ф) 0,2 


Первоначальное давление (Р) МПа 20 


Общая сжимаемость ( с, ) МПа" 10“ 


Рабочий дебит (4) м? / сут 100 


. Смоделированные данные ГДИС 
Модель течения в бесконечном пласте 


. Смоделированные данные ГДИС 
Модель течения в бесконечном пласте 


Изменение давления и его производная, МПа 
> 

Изменение давления и его производная, МПа 
2 


1 2 2 и 1 


10° 10" 10° 10 10° 10 10 10° 10 10? 


544... Ал... 


а) 51 точка данных в) 71 точка данных 


. Смоделированные данные ГДИС 
Модель течения в бесконечном пласте 


. Смоделированные данные ГДИС 


Изменение давления и его производная, МПа 
5. 
Изменение давления и его производная, МПа 


Модель течения в бесконечном пласте 
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Рисунок 1 — Смоделированные данные ГДИС 
и соответствие им корректно выбранной модели пласта 
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Таблица 3 — 95 %-е доверительные интервалы для проницаемости в случае кор- 


ректной модели 


а) 6) в) г) 
Число точек данных 51 61 71 81 
Оценка параметра 0,0484 0,0498 0,0498 0,0499 
52 2.74.10” 2,9010? 2.67.10” 2.69.10” 
п" 65510" 2,09.10° | 6,40.10* | 2,82.10“ 
ое="й 4,2410” 6,07 .10* 1,71.10° | 7,57.107° 
{7 6,5110 2.46.10" 1.31.10 8,70 .10? 
Доверительный интервал 2,71% 0,99 % 0,52 % 0,35 % 
Решение приемлемо приемлемо | приемлемо | приемлемо 


Для простоты приводятся результаты по оценкам только для одного параметра — 
проницаемости. Маргинальные распределения вероятностей изображены на рис. 2. Соответ- 
ствующие 95 %-е доверительные интервалы для проницаемости сведены в табл. 3. 


5000 т 


51 точка данных 
4500 скыниимины 61 точка данных | 
або 71 точка данных - 
81 точка данных 


Плотность вероятностей 


ь :. 
0.045 0.05 
Проницаемость, мкм 


2 


Рисунок 2 — Маргинальные плотности распределения вероятностей в случае 
корректно выбранной модели пласта 


По суги, только случаи 6), в) и г) содержат полезную информацию о 
проницаемости. Как следует из табл. 3, оценки проницаемости достаточно близки к 
истинному значению (А = 0,05 мкм’). Поэтому на рис. 2 все распределения вероятностей 
группируются вокруг этого значения. С увеличением числа данных появляется больше 
информации о проницаемости, и соответствующее ей СКО (с,) уменьшается. Размах 
распределений сужается, и нормальное распределение стремится принять форму дельта- 
функции Дирака. С позиций доверительных интервалов все случаи являются прием- 
лемыми, т.е. модель выбрана верно. 

Следующий случай соответствует некорректной с точки зрения данных модели. 
Данные ГДИС методом понижения уровня были получены с помощью модели непро- 
ницаемой внешней границы, на которые затем накладывались такие же, как и в пре- 


дыдущем случае, случайные ошибки. Истинные значения параметров А = 0,05 мкм”, 
$=10, С=0,2 м /МПа и Г =500 м. Для подгонки кривой изменения давления 


применялась модель течения в бесконечном пласте с тремя параметрами 
(А, БиС). Соответствие модели данным иллюстрируется на рис. 3. 
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Рисунок 3 — Смоделированные данные ГДИС и соответствие им некорректно 
выбранной модели пласта 


Маргинальные распределения вероятностей изображены на рис. 4. Соответствую- 
щие 95 %-е доверительные интервалы для проницаемости сведены в табл. 4. 


Таблица 4 — 95 %-е доверительные интервалы для проницаемости в случае некор- 


ректной модели 


а) 6) в) г) 
Число точек данных 51 61 71 81 
Оценка параметра 0,0484 0,0487 0,0362 0,019 
52 2.74.10? 4,77 10 6,59 .10° 1,98 -107' 
р' 1,55.10? 1,92 103 1,80 -10— 6,10.10° 
о = 4,24.107 9.14.10 118-10 т 
{7 6,51.10“ 3,0210“ 1.09.10? 110.103 
Доверительный 271% 1,24 % 6,00 % 11.49 % 
интервал 
Решение приемлемо приемлемо приемлемо | неприемлемо 
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Рисунок 4 — Маргинальные плотности распределения вероятностей в случае 
некорректно выбранной модели пласта 


В случае г) и д) модель течения в бесконечном пласте не соответствует данным. 
По сравнению со случаями а) и 6) средний квадрат ошибок 5’ для случаев в) и г) 
начинает возрастать, что приводит к увеличению размаха графиков плотностей 
распределения вероятностей. Если случаи а), 6) и в) еще приемлемы с позиций 
доверительных интервалов, то случай г) уже неприемлем, что говорит о некорректном 
выборе модели. В принципе доверительные интервалы можно использовать для 
принятия или отклонения выбранной модели. Независимо от того, корректно или нет 
выбрана модель, доверительные интервалы в конечном счете дают непротиворечивые 
результаты. Но нужно учитывать, что на практике при верификации модели 
доверительные интервалы следует определять для всех параметров. Кроме того, 
доверительные интервалы легко вычислять, поскольку вся необходимая информация 
содержится в результатах нелинейной регрессии, и нетрудно использовать для 
верификации модели, что было продемонстрировано выше. Однако сравнительный 
анализ на основе доверительных интервалов имеет два недостатка (практический и 
теоретический) с точки зрения дискриминантного анализа моделей. 

Во-первых, доверительные интервалы прямо пропорциональны дисперсии 
распределения вероятностей параметра, который в свою очередь представляет собой 
комбинацию среднего квадрата ошибок (оценочная дисперсия) 5’ и диагонального 


> —1 
элемента обратной матрицы Гессе й, . Например, в случае г) (табл. 4) оценка дисперсии 


составляет 6,59.10° МПа? , что намного больше истинной дисперсии 2,5.10” МПа?. 
Однако из-за того, что диагональный элемент обратной матрицы Гессе относительно 
мал (1,8.10*) вследствие большого числа точек данных (71 точка данных), соответ- 


ствующая дисперсия для проницаемости составила всего 1,18.10`°. В результате, 95 %- 
й доверительный интервал оказался приемлемым (6 %). Тем не менее, анализируя рис. 
3, можно предположить о существовании непроницаемой границы, т.к. нетрудно 
различить крутой наклон прямой на графике производной функции изменения 
давления. То есть доверительные интервалы могут находиться в допустимых пределах, 
даже если используется некорректная модель. 


330 «Искусственный интеллект» 42008 


Верификация модели пласта-коллектора по данным гидродинамических исследований 4 д 


Во-вторых, доверительные интервалы удобны при верификации моделей, но не 
подходят для их дискриминантного анализа. Другими словами, на основе довери- 
тельных интервалов можно определить, подходит ли модель или нет, но ничего нельзя 
сказать о том, какая из моделей лучше. Это связано с тем, что при вычислении 
доверительных интервалов не учитывается корреляция между параметрами. Однако в 
общем случае параметры пласта нелинейно связаны друг с другом, что необходимо 
принимать в расчет при верификации. Более того, уравнение (25) говорит о том, что 
размерность распределения вероятностей параметров совпадает с их количеством. То 
есть различные модели с различным числом параметров имеют различные размер- 
ности распределений вероятностей. Поэтому прямое сравнение соответствующих 
доверительных интервалов явно недостаточно. 
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С.В. Денисов, К.А. Сидельников 

Верифкащя модел! пласта-колектора за даними пдродинам!чних дослджень 

свердловин з використанням дов!рчих 1нтервал!в 

Представлена базова процедура нелийного регресйного анал1зу. Використаний байесовський 
висновок, при якому вся 1иформаця про параметри пласта виражаеться у терм!нах Ймов1рности, 
а невизначетсть, пов’язана з ощнками цих параметрив, може визначатися юльюсним чином. Досллджен! 
проблеми, як! виникають шд час використання дов1рчих 1нтервалв у межах байесовського висновку 
для дискримнантного анал1зу р1зних моделей-кандидатв. 


5.7. Реш5оу, К.А. У4аетйоу 

ТБе Базс ргосефигез оЁ попПпеаг гезгез$1оп аге ргеземе4. Вауезап ш{егепсе 1$ шнодисед. ш Вауезап 
шЕгепсе а] пЮгтайоп абойе Фе гезегуой рагате{егз 15 ехргеззе4 ш {ег оЁ ргоба цу, ап ипсематбу 
шуо[уе4 ш Фе рагате{ег езИтаез сап Бе ехргеззе4 диапщануе!у. СопНЯ4епсе пиегуа|5 аге депуе4 ш Фе 
Пате\уогК оЁ Вауезап шегепсе, ап фе ргоМетз шБегепИу шуо!уед ш @е аррйсаНоп оЁ сопН4епсе 
шбегуа!$ Юг тоде! Ч15сгиптайоп аге 415си$$е4. 
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